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ABSTRACT 

We compare the recently published velocity dispersions for seventeen Andromeda dwarf spheroidals 
with estimates of the MOND predictions, based on the luminosities of these dwarfs, with reasonable 
stellar M/L values, and no dark matter. We find that the two are consistent within the uncertainties. 
We further predict the velocity dispersions of another ten dwarfs for which only photometric data are 
currently available. 

Subject headings: dark matter — galaxies: kinematics and dynamics — Local Group 
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1. INTRODUCTION 

The Modified Newtonian Dynamics (MOND) 
(|Milgroml 11981 is a hypothesis that explains the mass 
discrepancies observed in galaxies without dark matter. 
It has recently been reviewed by iFamaev fc McGaughl 
(|2012f) . In MOND, the dynamics is conventional at 
large accelerations but modified at low accelerations 
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In the limit of very low 



accelerations, it is hypothesized that the acceleration is 
enhanced over the usual Newtonian value gw, tending 
towards a deep MOND limit a — » ^JgN^o for a "C aoQ 

The predictions of MOND have largely been confirmed 
in rotat ionally supported galaxies (for reviews and refer- 
ence see iSanders fc McGaughl2002t IFamaev fc McGaughl 
120121) . It has been less extensively tested in pressure sup- 
ported systems, in part because bright elliptical galaxies 
only reach the MOND regime at l arge radii where the 
data become sparse (see, however, Milgrom fc SandersI 
12001 IRichtler et aD l20ll IMilgroml l2012bT ). The n~ 
dictions of the theory become most acute in systems 
of low surface brightness, where the characteristic bary- 
onic surface density Sj, <C ao/27rG7, which corresponds to 
gN "C ao, and hence to a < an. 

The dwarf Spheroidal (dSph) satellites of the Milky 
Way and Andromeda provide a set of very low surface 
density, pressure supported stellar systems in which the 
predictions of MOND can be tested. This test was first 
applied to the cla s sical dwarfs of the Milky Way by 
iGerhard fc Spergell (fl992h . They found that while five of 
the seven dwarfs were in good agreement with MOND, 
two (Draco and Ursa Minor) required unpleasantly high 
mass-to- light ratios. T hese improved as the data im- 
proved ((Milgroml (1995) , but persist in being marginal 
cases when analyzed with only a global velocity disper- 

1 This formulation occurs naturally if dynamics is required to 
be invariant under space-time scaling, (t, r) — > A(t,r). In self- 
gravitating systems, the Newtonian force per unit mass </jv ~ 
MG/r 2 scales as A -2 , while accelerations, a, scale as A -1 (M is 
the observed luminous mass) . We can then not have a = , which 
is not scale invariant. Instead, a relation of the form a ~ v / ab~9jv 
is implied. For any isolated mass we must have asymptotically far 
i = Q^/ a 09N> with a universal constant q; ao is normalized so that 
g = l {Milgrom 2 0~T2aT) 



sion ((McGaugh & Wolf 2010]). One must always bear in 
mind the applicability of the necessary assumptions; in 
Ursa Minor for exampl e there is eviden ce for multiple 
dynamical components ((Pace et alll2012[ ) which, if true, 
would considerably complicate the analysis. However, 
the more basic issue was whether we should be more im- 
pressed by the five cases that are in good agreement with 
MOND or the two that disagree. 

The data for the c lassical Milky Way d warfs are consid- 
erably better now. I Walker et al.l (|2007t ) have measured 
the velocity dispersion profile as a function of radius. 
This enables detailed analysis of their ma ss profi l es, not 
just a single global ch aracteristic value. lAngusI (2008) 
and ISerra et ail (|2010f ) found these to be in good ac- 
cord with MOND provided that the orbital anisotropy 
varies with radius , much as found in conventional anal- 
yses ((Evans et al.l 120091 ) . At this level of detail, one 
must also be cog nizant of tidally liberated interlopers 
(|Serra et aLll2010l) , and differing velocity dispersion pro- 
files for diff erent stellar populations (e.g..,Battagli aet al.1 

EonSdlp. 

In addition to improved data, the available sample 
has grown rapidly with the discovery of many new 
dwarf satellites around both the Milky Way and An- 
dr omeda. The so-cal l ed ult rafaint dwarfs are discussed 
bv iMcGaugh fc WolJ ((20101) . who find that these tiny 
galaxies have velocity dispersions larger than expected 
for isolated systems of their mass in MOND. However, 
these are not isolated systems. Tidal effects due to the 
Milky Way cannot be ignored, obviating the usual as- 
sumption of dynamical equilibrium. It will be interest- 
ing to see if any of the ul trafaints are in the proce ss of 
dissolution anticipated bv lBrada fc Milgroml ((20001 ). 

An ear ly test of MOND with the dwarfs of M31 was 
made bv ICote et al.l (|1999T) . who analyzed And II. For 
a measured velocity dispersion of a = 9.3 kms -1 , they 
found good agreement with MOND with no dark matter: 
M*/Lv = 3.2 M /L Q instead of M/L v = 21 M Q /L in 
the Newtonian analysis. Differing velocity dispersions 
have since been publis hed, and the estimated luminosity 
of Ar id II has tripled (|Kalirai et a,l.ll2010t iTollerud et all 
2012), lowering M * /L accordingly. We extend the 
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MOND analysis to the many more recently measured 
dwarfs (and updated data for And II). 

There are at present seventeen Andromeda dwarfs with 
published velocity dispersion measurements that can be 
employed for this test. There are a further ten that do 
not yet have published velocity dispersions but for which 
there exists sufficient photometric information that a pre- 
diction can be made. 

These dwarfs all exhibit large mass discrepancies, re- 
quiring most of their mass to be dark in Newtonian 
dynamics. On the other hand , the recent realization 
(jlbata et all 120131 : iConn et al.l 120131 ) that many (15) 
dwarfs (including twelve that we analyze here) are part 
of a large, kinematically coherent disc of satellites, is at 
odds with the notion that they are the primordial oc- 
cupants of the sub-halos found in ACDM simulations: 
their phase space distribution is quite different. Sub- 
halos comprise a dynamically hot, quasi-spheircal popu- 
lation rather than a thin, rotating disk. Arguments to 
the same effect have been made in connection with a 
similar disk of satellites around the Milky Way (iKroupa 
120121: iPawlowski et al.ll2012l ). iDabringhausen fc Kroupa 
(2013) suggest that these orbital properties are the re- 
sult of formation as tidal dwarfs. Tidal dwarfs should 
not exhibit the observed mass discrepancies in ACDM 
because of the difference in pha se space between dar k 
matter and baryonic tidal debris (jBournaud et al.ll2007D. 
but this situation can arise naturally in MO ND ( M ilgroml 
[20071 : IGentile et al.ll2007t iTiret fc CombesH2008[ ). 

In MOND we predict at the outset that these dwarfs 
should show large mass discrepancies, as observed. This 
is based on the strong correlation between low acceler- 
ation and high mass discrepancy predicted by MOND 
(Milgrom 1983). These low surface brightness systems 
all have very low internal accelerations. The accelera- 
tion at the half mass radius is estimated to be between 
O.Olao and 0.15ao for all dwarfs considered here, with 
only three exceeding 0.1a . 

In $2] we describe the MOND estimators we use for 
the bulk velocity dispersions. In Sj3]we discuss the various 
data we use. In $4j we describe our results. These include 
a comparison of the observed velocity dispersions with 
those predicted by MOND for seventeen dwarfs f ^4.1[) . 
For another ten dwarfs, we predict the velocity dispersion 
in advance of measurement (i j4.2l) . We conclude in $5] 

2. MOND ESTIMATES OF THE VELOCITY 
DISPERSION 

Our basic estimate of the predicted MOND line-of- 
sight velocity dispersion is based on the deep-MOND 
virial relation 

MGa Q = \{v 2 ) 2 . (1) 

Here, (v 2 ) is the 3-D mass weighted mean-square veloc- 
ity, and M is the total mass of the system. This relation 
holds exactly in t he deep-MOND limit of the m odified 
Poisson equation of Bckcnste in fc Milgroml (|1984[ ) for iso- 
lated, self-gravitating, stea dy-state (virialized) s ystems 
made of test particlefl (see iMilgroml 11994 119971) . It is 

2 Namely, constituents whose masses are all much smaller than 
the total mass of the system. There is a simple g eneralization of 
this relation to systems made of arbitrary masses HMilgromlll9971 . 
|20TC|) . 



also exact in the QUM OND formulatio n of MOND, un- 
der the same conditions (Milgrom 2010). More generally, 
for any isolated, deep-MOND, self-gravitating system, a 
relation of the form Q 2 MGao — V 4 should hold for any 
MOND th eory (as it follow s from the basic assumptions 
of MOND: lMilgrom1l2012a| ). where V is some measure of 
the mean 3-D velocity, and Q ~ 1 may depend on the 
exact theory, and on details of the system. 

The available data on the M31 dwarfs give only some 
measure the line-of-sight velocity dispersion, a, which 
we want to compare with (v 2 ) 1 ^ 2 deduced using eq. ([T]). 
For the comparison to be meaningful we have to assume, 
in default of better knowledge, that the M31 dwarfs 
are globally isotropic, namely that a would be indepen- 
dent of the line of sight. This would not be quite cor- 
rect if rotation is present, or if the system is otherwise 
globally anisotropic. Furthermore, we assume that the 
quoted a represents the mass-weighted dispersion over 
the whole system, which may not be quite correct in gen- 
eral. Specifically, in the case of the M31 dwarfs a is based 
only on measured velocities of red giants, and the spatial 
coverage of each system is not always complete. It is not 
clear to what extent the necessary assumption that the 
observed velocity dispersion of a particular stellar popu- 
lation along our line of sight is a valid approximation of 
the desired 3D mass weighted velocity dispersion. 

Making all those assumptions permits us to represent 
(v 2 ) by 3cr 2 , and thus to write the line-of-sight dispersion 
predicted by MOND as 

Visa « (^MGa ) 1/4 . (2) 

We use the subscript iso to specify the case of an isolated 
object not influenced by external fields. 

The assumption of isolation is clearly violated for some 
of the M31 dwarfs. Isolation here means that the sys- 
tem is not subject to a n appreciable external-field effect 
(EFE) (|Milgromlll983l) . Let j m represent accelerations 
within the system (here, the dwarf) with respect to its 
center of mass, and g ex the acceleratior0 with which the 
system falls as a whole — here the acceleration produced 
by M31 at the position of the dwarf. Then, relation ([2]) 
applies only if gi n 3> g ex . Some of the M31 dwarfs sat- 
isfy this condition to a reasonable approximation. In the 
opposite case, gi n <C g ex , the EFE is dominant. In such 
cases, the above mentioned MOND theories predict that 
the internal dynamics of the system is quasi-Newtonian: 
the Newtonian virial relation is satisfied with two cor- 
rections: i) instead of Newton's constant G one has to 
use a larger effective value G e ff = G/n(g ex /ao), where [i 
is the MOND interpolating function, and ii) the dynam- 
ics is not rotationally invariant, since the external field 
introduces a preferred direction. We thus expect some 
induced anisotropy. In the case of EFE dominance, our 
estimator becomes 



3 Here gj n is taken as the acceleration on a star at the half- 
mass radius due to the mass of the dwarf itself, while g ex is the 
acceleration on the center-of-mass of the dwarf due to M31. These 
are the actual accelerations, not just the Newtonian expectation 
<7iV f° r the baryons alone. 
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where rw2 is the 3D half light radius of the dwarf. 

It so happens, however, that the M31 dwarfs that are 
not relatively free of the EFE, are also not of the other 
extreme type: they have gi n ss g ex . It is then rather 
difficult to make an estimate of the MOND prediction for 
a. For these cases we simply give the predictions MOND 
would make for a, for both the isolated case according 
to eq. ([2]), and the EFE dominance case of eq. ([3]). The 
predicted value of a for a given mass should always be 
smaller when the EFE is important, provided the system 
remains in dynamical equilibrium. Systems subject to 
very strong external fields may suffer tidal stirring and 
disruption. 

For each dwarf we estimate g ex = V^ 3 1 /Dm31, where 
V M 3i ~ 230 kms" 1 JCh emm etaLl 2009) is an estimate 
of the asymptotic rotational speed of M310, and £>M3i 
is the distance from the dwarf to the center of M31. We 
estimate gi n = 3<r 2 /ri/ 2 , where a is calculated from eq. 
((2|) . These are used to determine the importance of the 
EFE, and g ex is also used in calculating G e //. In the 
limit of very small accelerations that is satisfied here, 
deff ~ Gao/g ex . 

3. THE DATA 

Spectra of individual stars in dwarf spheroidal galax- 
ies are used to construct the velocity dispersions of their 
stellar population. At the distance of M31 (~ 780 kpc), 
only red giants are amenable to observation. It is pos- 
sible for different stellar populations to have different 
velo city dispersions; this situation is o bserved in For- 
nax (iBattaglia et al.ll2006l ) and Sculptor (Battagli a et al.l 
2008). In order to perform the analysis here, we make 
the usual assumptions that the observed, line-of-sight ve- 
locity dispersion of the red giants is representative of the 
full mass-weighted 3D velocity dispersion and that the 
system is in dynamical equilibrium. 

The assumption of dynamical equilibrium is not as safe 
in MOND as it is in the context of dark matter. In the 
latter case, the stellar extent of dwarf spheroidals is usu- 
ally well within their inferred tidal radii: the stars are 
safely ensconced within their dark-matter halos. Tidal 
effects are relatively stronger in MONE0, and there is no 
protective cocoon of dark matter. To complicate matters, 
a dwarf with an eccentric orbit that brings it close to its 
host may suffer non-adiabat ic perturbations at pericen- 
tcr (|Brada fc Milgroml[2000l ) that may persist even as it 
spends most of its time near apocentei@. We do not know 
the orbits of the dwarfs of Andromeda, only their present 
distances from M31. At their present locations, only And 
IX appears to be in imminent danger from tides as it has 
-^/ r i/2 < -^M3i/-Dm3i (where M is the stellar mass of 
the dwarf and Mm31 the baryonic mass of the host). The 

4 The predicted EFE velocity dispersion <T e j- e depends in- 
versely on Vm31, so would be 10% lower if Vmh were 10% high 
(253 kms -1 ). The approximation is necessarily crude, as for some 
of the more distant dwarfs, the acceleration due to other objects 
may begin to play a role (e.g., the effect of M33 upon And XXII, 
where the opposing vectors of the larger galaxies may tend to can- 
cel one another and reduce the EFE). 

5 Ursa Minor, Draco, and nearly all of the ultrafaint dwarfs 
have photometric tidal ra dii comparable to their MOND tidal radii 
(|McGaugh fc Wolfll201Ch . 

6 Perhaps th e tidal debris recently observed around Carina by 
Battaglia et al. (2012) is an example of this effect. 



rest of the dwarfs of M31 are sufficiently far from their 
host that they should not suffer tidal disruption unless 
they are on rather eccentric orbits. 

The luminosities, half light radii, velocity dispersions, 
and M31-centric distances employed here come from a 
variety of sources. The most exten sive kinematic sur- 
vey to date is the SPLASH survey (jKalirai et al.l 120101: 
iTollerud et al.l [20121) . In addition to providing a large 
and fairly homogeneous set of data, dwarfs for which 
there are multiple analyses illustrate how the deduced ve- 
locity dispersion can change with further observation or 
simply different analysis. Usually the results are consis- 
tent, but in some cases the velocity dispersion changes by 
more than the formal uncertainties would predict (com- 
pare, for example, published estimates of the velocity 
dispersions of And II, And III, and And VII ). 

Further data are pro vided bv lCollins etaLl (|20inl . l2f)Tlh 
and iChapman et al.l ()2012l ). In addition to providing 
measurements of dwarfs not yet covered by SPLASH, 
these authors also provide completely independent mea- 
surements of some of the same dwarfs, providing another 
consistency check. There are at present four cases of 
dwarfs with overlapping, independent data. In the case 
of And V, the measured velocity dispersions are con- 
sistent. In And IX, the velocity dispersions differ sub- 



stantially: g = 10.9 kms vs. 4.5 kms 
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case of And XIII there is also a difference in the mea- 
sured v alue, but the formal sig nificance of this difference 
is less. IChapman et al.l (|2012| ) estimate an upper limit 
for the velocity dispersion of An d XXII that is consis- 
tent with the value measured by ITollerud et al.l (|2012|L 
Differences between independent studies, when present, 
may be attributable to different observational techniques 
(e.g., fibers vs. slits), but do not strike us as beyond 
the realm of happenstance given our experience with the 
dwarf satellites of the Milky Way. 

Table [T] compiles the relevant observational data. Col- 
umn 1 gives the name of the dwarf. There are gaps in the 
numerical sequence when data do not exist or the rele- 
vant object turns out to be a backgroun d galaxy rather 
than a satellite of M31 (e.g., And IV: iFerguson et al.1 
2000) . The second column is the luminosity of each dwarf 
and the third its 3D half light radiusQ The fourth and 
fifth columns give different measurements of the observed 
velocity dispersion. Examination of these two columns 
permit one to judge the consistency of repeat measure- 
ments. The last column gives the references from which 
these data were taken. The remaining columns present 
the predictions of MOND as discussed in the next sec- 

7 The 3D half light radius has typically been deprojected from 
the observed half light radius (r-1/2 = 4i? e /3) assuming spherical 
syste ms with velocity dispersions independent of radius IWolf et al.l 
1201011 . The latter assu mption seems reason able given the example 
of the classical dwarfs (Walker et al. 2007), but is not guaranteed 
to hold. The assumption of sphericity is a more subtle issue. The 
quoted half light radii are typically major axis values, appropriate 
for oblate systems. For prolate or triaxial system, the geometric 
mean of the axis ratios might be a better proxy for r 1 /2- Since 
the half-mass radius only enters the MOND calculation for the 
EFE case, and then only as the square root, we simply adopt the 
published major axis values for consistency with previous work. 
In the most flattened case (And XXII with e = 0.56), adopting 
the smaller geometric mean for the radius would result in a ~ 
7% increase in cr e f e . At this level, anisotropy is also likely to be 
important, and cause the observed line-of-sight velocity dispersion 
to deviate from the predicted mass- weighted 3D velocity dispersion. 
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Table 1 

Predicted and Measured Velocity Dispersions 



Dwarf 


L v 


r l/2 


°of>s,l 


a obs,2 


G iso 




9 1 n 


Ref. 




10 5 L 


pc 




km s 






ao 




And I 


45. 


832 


10.2 ± 1.9 


10.6 ± 1.1 






0.082 


1,2 


And II a 


93. 


1660 


7.3 ±0.8 


10.0 ± 1.7 


n.o±H 




0.059 


2,3 


And III 


10. 


525 


9.3 ± 1.4 


4.7 ± 1.8 


6.3±|| 


5-4+?' 3 6 


0.062 


1,2 


And V 


5.9 


292 


10.5 ± 1.1 


ii.5±|| 


^'^ — 0.9 


7.0+ti 

— 2.1} 


0.085 


1.1 


And VI 


34. 


440 


— 2.4 




8.6+H 

— 1.4 


0.14 


4 


And VII 


178. 


977 


13.0 ± 1.0 


9.7 ± 1.6 


13.0+2? 

— 2.1 




0.14 


1,2 


And IX b 


1.5 


552 


10.9 ± 2.0 


— 3.4 


3 n+0.7 
— 0.6 


1-5+21 

—0.4 


0.023 


1,5 


And X 


0.76 


309 


6.4 ± 1.4 


3.9 ±1.2 


o q+0.6 
0.5 


9 Q + 1.0 

z -°-0.7 


0.029 


1,2 


And XI 


0.49 


145 


< 4.6 




3 n+0-6 
a,u -0.5 


q 9 + 1.3 
J ' z -0.9 


0.049 


5 


And XII 


0.31 


289 


2 6+ 51 


q 7 +8.9 

y -'-4.5 


2 6+ - 5 


1 7+0.7 

1 -' -0.5 


0.020 


5 


And XIII 


0.41 


203 


5.8 ± 2.0 


9 o+0.5 


2 5+ 10 


0.032 


1,5 


And XIV 


2.1 


537 


5.3 ± 1.0 


5.4 ± 1.3 


ft -°— 0.7 


3 g+1.5 


0.027 


1,2 


And XV 


7.1 


355 


4.0 ± 1.4 




5 - 8 -0.9 


6 1+ 2 ' 5 
°' i -1.8 


0.076 


1 


And XVI 


4.1 


178 


3.8 ± 2.9 




5 n+ 10 




0.12 


1 


And XVIII 


6.3 


417 


9.7± 2.3 




5 6+ 11 
°-°-0.9 




0.061 


1 


And XXI 


1.6 


1023 


7.2 ± 5.5 




k o+l-O 
°- z -0.8 


q 7+1.5 
-'-l.l 


0.021 


1 


And XXIP 


0.35 


340 


3-5±|;l 


< 6.0 


9 7 +0.5 
z -'-0.4 


9 q+1.0 
^•°-0.7 


0.018 


1,6 


References. 


— The luminosities, half light radii, and velocity dispersions are 


taken 



from 1. Tollcrud ct al. (2012), 2. Kalirai ct al. (2010),' 3. Ho ct al. (2012), 4. Collins ct al. 
l|20"TlT ). 5.|Collins ct al. (2013), and 6. IChaprnan et all pOilh . 



Note. — The velocity dispersion predicted by MOND in the isolated and EFE limits 
is given for M*/L = 2±j M /L . The correct regime depends on M*/L; one should 
generally take the lower of the two predicted velocity dispersions, so there EFE column is 
omitted when it exceeds the isolated prediction for the range of considered mass-to-light 
ratios. The internal acceleration at the 3D half mass radius is estimated in units of ao 
assuming the isolated case with M„/L = 2 M /L0. 

a And II is unusual for dSph galaxies in having measured rotation (Ho et al. 2012). This is 

not included in the older cr obg l , but is included in <r oi)s 2 through cr obs 2 = \J ® 2 + \ v 2 - 

b The predictions for And IX are very uncertain as this dwarf is subject to a strong external 
field from M31. This makes the assumptions of isotropy and dynam ical equilibrium suspect . 
c We consider only the EFE due to M31, where we use ri/ 2 from[Chapman ct al. (2012). 
The larger radius given by Tollcrud ct al. (2012) would give cr e f e smaller by a factor 0.8. 
IChaprnan et all) (2012) argue that And XXII might be a satellite of M33 rather than M31. 
The accelerations due to M31 and M33 acting on And XXII are quite uncertain, but appear 
to be approximately equal in magnitude yet opposite in direction. Vector cancelation of 
these external fields might have the net effect of reducing the importance of the EFE. 



tion. 

4. RESULTS 

We use the formalism described in § [5] to calculate 
the predicted velocity dispersion for each of the dwarf 
spheroidals of M31, using the photometric datc@ in Ta- 
ble [TJ The results are given in Table [1] and shown in 
Fig. [TJ together with the different measured velocity dis- 
persions. For isolated objects in the deep MOND regime, 
only the stellar mass is required. For objects where the 
external field dominates, the half-mass radius is also nec- 
essary, appearing as usual in the conventional virial the- 
orem with only the effective value of 67 altered (§ [2]). 
We adopt the specific valu^3 ao = 1-21 x 10~ 10 ms~ 2 

8 Where published photometric data differ, we adopt the non- 
SPLASH values as a hedge against the systematic uncertainties 
that might afflict any one sample. Affected objects are And V, IX, 
XIII, and XXII. The differences are minor for all but And XXII, 
where r]y 2 differs by a factor of ~ 2 and the luminosity differs by 
~ 50%. 

9 The value of ao obtained by Bcgcman ct al. (1991) remains 
consistent with more recent estimates (McGaugh 2011). 



(|Begeman et al.|["l991|) to perform this calculation so that 
(cr/kms" 1 ) = (M/1264MQ) 1 / 4 . 

The luminosity is typically more accurately measured 
than the velocity dispersion, so we use this to predict the 
velocity dispersion with MOND for a range of assumed 
mass-to-light ratios [M — {M*/L)L). We compute a for 
M^jLy = 1, 2, and 4 M q /Lq to show the sensitivity of 
the result to the assumed mass-to-light ratio. There is no 
guarantee that the stellar population of any given dwarf 
is within these limits, but it does span a plausible range. 
It also gives an impression of the range in the prediction 
and its dependence on the stellar mass-to-light ratio. 

The first calculation assumes that each dwarf is iso- 
lated. The isolated MOND prediction <7i so is tabulated 
in column 6 of Table [JJ where it can be compared directly 
to the observations. The uncertainty in the prediction 
only represents the above range of mass-to-light ratios. 
It does not represent the uncertainty in any observed 
quantity nor does it include any of the potential system- 
atica discussed above. Here the luminosity is the only 
observational input, for which the uncertainty should be 
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Figure 1. The predicted and observed velocity dispersions of the dSph satellites of M31. Each bin along the x-axis represents the 
corresponding dwarf: And I is 1, And II is 2, and so on. Round symbols are the predictions of MOND. Large symbols assume Mt/Ly = 
2 Mq/Lq. Small symbols represent M*/Ly = 1 Mq/L© at the lower end of the predicted range and M*/Ly = 4 Mq/Lq at the upper 
end. The calculation in the isolated case is shown as a filled green circle; that when the EFE dominates is shown as an open circle. The 
appropriate case is shown for each mass-to-light ratio based on which p redicts the lower vel ocity dispersion, as disc ussed in the tex t . Dat a 
from the SPLASH survey are s hown as squares: open squ ares are from iKalirai ct al. (2010); filled squares are from Tollcrud ct al. (2012). 
Triangles represent the data of ICollins et al.l (120101 . 120111 ). In the case of And II, the earl y measurement of [Dote et aTT " iIT999l ) is shown as 
an X , while the recent detection of rotation i n this case is included in the measurement of lHo et al.l (120121 . star) . The recent upper limit in 
the case of And XXII (Chapman ct al. 2012) is shown as the downward pointing triangle. 



less than that in M*/L. We simply tabulate <Ti SO for 
M^/Ly = 2 M Q /L and give the range for M*/Ly = 1 
to 4 Mq/L© as the uncertainty. 

We confirm that the dwarfs are indeed in the MOND 
regime by estimating the internal acceleration at the half 
light radius, again assuming M*/Ly = 2 Mq/L©. This 
is given in column 8 of Table [T] in units of ao . All of 
the dwarfs considered here have gi n < 0.15ao, and of- 
ten considerably less, as expected for systems of low sur- 
face brightnes£3. This places them deep in the MOND 
regime where the details of the interpolation function are 
irrelevant. 

While all dwarfs are in the deep MOND limit, it is less 
obvious that they qualify as isolated systems. We check 
whether the assumption of isolation holds by estimating 

10 The connection between surfac e density and acceleration is 
axiomatic in MOND (Milgr oml IT983t ). Whether it holds in ACDM 
depends on the specifics of galaxy formation theory. 



the external acceleration g ex upon the dwarf due to M31, 
and comparing that to the internal acceleration g^ n . As 
is the case for the satellites of the Milky Way, g ex is often 
comparable to gi n . In many cases, the two are the same 
within the errors. 

We therefore also report the velocity dispersion for the 
EFE dominated case in column 7 of Table [TJ In most 
cases, the values for M^/Ly = 2 M©/L© are not much 
different than in the isolated case. In the deep-MOND 
limit the ratio of the dispersions calculated for the two 
cases is <7 e /e/ciso = y^gm/gex- Consequently, the case 
that reports the lower velocity dispersion is always the 
formally correct one. We therefore report cr e / e in Ta- 
ble Q] only if it is predicted to be less than that in the 
isolated case for at least one choice of mass-to-light ra- 
tio. The range in predicted velocity dispersion spanned 
by the same range of 1 < M^jLy < 4 is larger than in the 
isolated MOND case, owing to the Newtonian a cx M 1 / 2 



() 



scaling] rather than the MOND a oc M 1 / 4 scaling. Here 
too we have not included the uncertainties due to dis- 
tance and luminosity uncertainties^ 

Any one dwarf may transition from the isolated to EFE 
regime for different choices of mass-to-light ratio. This 
is denoted in Fig. Q] by marking each point as a filled 
or open symbol for isolated or EFE cases, respectively. 
One should bear in mind that the assumption of isotropy 
is immediately suspect in cases dominated by the EFE, 
which can induce anisotropy. The EFE predictions are 
thus not as robust as those for isolated dwarfs, and we 
tabulate for reference <7i SO in all cases, even when a dwarf 
appears to be entirely in the EFE regime. 

4.1. Individual Cases 

Our predictions are most secure for isolated dwarfs. 
Based on the calculation described above, dwarfs that 
appear to be isolated for all choices of mass-to-light ratio 
are And II, And VI, And VII, And XVI, and And XVIII. 
And V is also likely to be in this category, though for 
M*/Lv = 1 Mq/Lq it becomes margina lly EFE domi- 
nated if we adopt the photo metric data o f Tollerud et al.l 
(|2012D rather than those of iCollins etail ([201 ID . Dwarfs 
that appear to be dominated by the EFE for all choices 
of mass-to-light ratio are And IX, And X, And XII, 
And XIV, and And XXI. The others span both regimes. 
This designation only reflects the assumed variation in 
M*/Ly; in many cases the observed internal acceleration 
3 a obs / r i/2 is difficult to distinguish from the estimated 
external acceleration V^ 31 / Dm3i given the observational 
uncertainties. 

And V is the worst case among the present data. It 
has a predicted velocity dispersion of <7j SO = 5.5 kms -1 
for M*/Lv — 2 M0/L0 compared to 10.5 ± 1.1 from 
ITollerud et all (l20ll and 11.5±^ from ICollins et al.1 
(|2011[) . In terms of mass-to-light ratios, the observed 
velocity dispersions imply M^jLy sa 10 Mq/Lq, im- 
plausibly high for a stellar population. The important 
question, however, is how significant the differ ence is. 
The stated uncertainty of ITollerud et al.l (12012D is con- 
siderably smaller than that of ICollins et al.l ([2011). The 
velocity dispersion we compute for M*/L y = 2 M^/Lq 
is 4.5 a smaller than the measurement of ITollerud et al.l 
(120121) but o n ly 1.4 a smaller than the measurement of 
Collins et al. (2011). Given the nature of the observa- 
tional data, and the way in which the velocity dispersion 
can change from observation to observation, it is unclear 

11 As in the isolated case, the uncertainty only represents the 
range in plausible mass-to-light ratios and not any observational 
or systematic errors. There is more opportunity for the latter in 
the EFE dominated case since we now need to know the rotation 
curve of M31 to large distances and the M31-centric distance of 
each dwarf. 

12 Uncertainties in the distance D from us to the dwarf are gen- 
erally not so large, since the dwarfs are roughly at the distance 
of M31. Since D enters the predicted values of <r iso in a square 
root, uncertainties in cr; so due to those in D are expected to be 
small compared with those due to uncertainties in M/L. How- 
ever, distance uncertainties are much more important for the pre- 
dicted <x e y e . Distance still enters cr e f e in a square root [through 
(M/r 1/ / 2 ) 1//2 ]i but, more importantly, errors in D can be greatly 
magnified as errors in there M31-centric distance Dm31- This en- 
ters <J e f e in a square root [(r e f e scales as (DM31D) 1 / 2 ] , besides 
entering the isolation criterion itself. For some dwarfs these uncer- 
tainties are quite significant, e.g., for And IX. 



how much weight to place on the stated uncertainties. 

The predicted velocity dispersion is too low in both 
And XIII and And XVIII, though here the significance 
is less (< 2a). In the case of And XIII, there are two 
independent observations, both higher than predicted 
albeit by a modest amount (~ 1.6a). There is only 
one observed value of the velocity dispersion in And 
XVIII at present, so there is no independent indica- 
tor of the accuracy of the uncertainty. Taken at face 
value, the predicted velocity dispersion is 1.8a too low 
for M*/Lv = 2 Mq/Lq. These do not seem like par- 
ticularly worrisome discrepancies given the nature of the 
data and the necessary assumptions in the prediction. 

And IX is an interesting case. The predicted dis- 
persion is low but consistent with the measurement of 
Colli nset al.l (12010J). It is substantially lower than that 
of ITollerud et all ([20121) . The data differ, so the inter- 
pretation depends on which measurement happens to 
come closer to reality. However, there is a further com- 
plication in the case of And IX. Of all the dwarfs dis- 
cussed here, this object is most clearly subject to strong 
tides from M31 based on the criteriE^ discussed by 
IBrada feTM ilgrom (2000]). This can result in the induced 
anisotropy mentioned in §[5J Depending on its orbit, And 
IX may also be subject to tidal disruption and dissolu- 
tion along the li nes inferred for the ultr afaint dwarfs of 
the Milky Way (IMcGaueh fc Woial20lol) . Consequently, 
the assumptions that underpin our analysis may be vio- 
lated, making our prediction for And IX the least secure 
of all the Andromeda dwarfs. 

Another interesting case is And II. ICote et al.l (1999) 
made an early velocity dispersion measurement of 
9.3^2'g kms -1 that is in good agreement with the 
MOND predict ion. This ap p eared to change with the 
observations of IKalirai et all (|2010f) . who found a lower 
velocity dispersion with a much smaller uncertainty: 
7.3 ± 0.8 kms" 1 . The MOND prediction for M*/L v = 
2 Mq/Lq is 4.6a higher than this sub sequently measured 
value. However, IKalirai et all (|2010f) were only able to 
measure one side of the dwarf, leaving open the possibil- 
ity that the difference in spatial coverage was the cause of 
the apparent difference i n mea sured velocity dispersions. 
Very recently, IHo et all (|2012l ) have reported substantial 
rotation in And II, with v = 10.9±2.4 kms -1 in addition 
to a dispersion of a = 7.8 ±1.1 kms -1 . This makes And 
II unique among dwarf spheroidals in having substantial 
rotational as well as pressure support. We incorporate 
both rotation and dispersion through a = \/ of os + I3v 2 
with j3 = 1/3 since ordered rotation enters the 3D rms 
velocity of eq. Q] in quadrature. This prescription is not 
exact, as the plane of maximum rotation is unknown, and 
the maximum rotation speed is not the mass weighted 
mean. These two effects tend to counteract one an- 
other, so this is probably the best available estimator 
(Wolf et al.ll2010h . Certainly it is an improvement on es- 
timates made in ignorance of rotation. Taking rotation 
into account, the data are in good agreement with the 
prediction of MOND. 

Indeed, the bulk of the sample is in good agreement 

13 And IX is currently within the MOND-prescribed Roche limit 
of M31. How it responds to this strong external perturbation de- 
pends on its orbit: how close it comes to M31 and how long it 
spends so close. 
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Table 2 

Predicted Velocity Dispersions 



Dwarf 


Ly 


r l/2 




a efe 


fliri 




10 5 Lq 


pc 


kms 


— i 


ao 


And XVII 


2.6 


381 


4 5+0-9 


z -°-0.7 


0.043 


And XIX 


1.1 


2244 


°- u -0.8 


^•"-0.8 


0.009 


And XX 


0.28 


165 




o 1+0.9 

^-o.e 


0.033 


And XXIII 


10. 


1372 




44+1.8 


0.024 


And XXIV 


0.94 


489 


o 5+O.7 


^•°-0.8 


0.020 


And XXV 


6.5 


945 


5 7+ 1 - 1 


■ J -°-1.0 


0.027 


And XXVI 


0.59 


296 


o 1+0.6 
d - i -0.5 


2 0+ - 8 
z - u -0.6 


0.027 


And XXVII 


1.2 


579 


o 7 +0.7 
°-'-0.6 


1 8+ - 7 


0.020 


And XXVIII 


2.1 


284 


4 o+O.S 




0.053 


And XXIX 


1.8 


481 


4 1+0-8 


3.8±i;S 


0.028 



Note. — Dwarfs listed here currently lack published 
kinematics. The velocity dispersion predicte d by MOND is 
based on the photometric data in Table 3 of McConnachie 
p012T ). We have taken the 3D r 1/2 to be 4/3 of the 2D half 
light radius reported by McConnachie (2012) for consistency 
with Table [T] Columns are the same as in Table [T] 

with MOND. The observed and predicted velocity dis- 
persion coincide for And I, And II, And III, And VI, 
And VII, And X, And XI, And XII, And XIV, And XV, 
And XVI, And XXI, and And XXII. For the cases in the 
isolated, deep MOND limit where our prediction is most 
secure, agreement is frequently bang on: compare 0~ o b s 
with Uiso for objects lacking a e f e in Table [1] And XVIII 
is a bit off, albeit by < 2a as discussed above. Only And 
V disagrees with MOND by as much as 4tr, and only then 
for one of the two measurements for that obi ect . The 
veloci ty dispersion reported for And II by iKalirai et al.l 
(2010) appeared to be problematic, but the detec t ion o f 
some rotational support in this dwarf bv lHo et al.l (|2012l ) 
brings it into agreement with the prediction of MOND. 

4.2. Further Predictions 

iMcConnachid (|2012| ) provides photometric information 
for ten dwarfs for which no velocity dispersion is cur- 
rently available. We can use this information to predict 
their velocity dispersions. MOND predictions for And 
XVII, And XIX, And XX, And XXIII, And XXIV, And 
XXV, And XXVI, And XXVII, And XXVIII, and And 
XXIX are given in Table [2] for the range 1 < M*/Xy < 
4M Q /L Q . 

These objects are deep in the MOND regime. Being 
quite faint, they are predicted to have low velocity dis- 
persions: a < 5km s" 1 for M*/L v = 2 M Q /L Q . The 
a priori predictions given in Table [2] can be tested as 
kinematic measurements become available. 

5. CONCLUSIONS 

MOND appears to be in good agreement with the ob- 
served velocity dispersions of the dwarf spheroidals of 
M31. Simply looking at Fig.Q] one might easily mistake 
the MOND predictions for a set of observational data. 

There is no a priori guarantee that MOND will work 
in these systems. In the context of the dark matter 
paradigm, we might observe dwarfs that reside in the 
~ 30 km s -1 sub-halos commonly se en in ACDM sim- 
ulations ()Bovlan-Kolchin et al.l l20lll) . or tidally formed 
dwarfs devoid of dark matter (jBournaud et al.l 12007ft 



with equilibrium Newtonian velocity dispersions < 
2 kms -1 . These are plausible scenarios with dark mat- 
ter that would result in velocity dispersions that MOND 
could not explain. 

Given the photometric data, the velocity dispersion 
follows in MOND. The value of an. is fixed to that es- 
tablished long ago fo r much brighter rotating galaxies 
(|Begeman et al.Nl99lD . The only computational flexibil- 
ity is in the stellar mass-to-light ratio, which is unavoid- 
able in any theory For a modest (± factor of two) range 
of M*/L motivated by our expectations for stellar popu- 
lations, the predictions of MOND are immediately in the 
right vicinity. 

As the velocity dispersion data continue to im- 
prove, it may be possible to infer M*/L for individ- 
ual dwarfs. We can then ask whether the MOND- 
inferred mass-to-light ratios correlate with indicators 
of population age observed in color-magnitude dia- 
grams. This is already claimed to b e the case in the 
dSph population of the Milky Way (|Hernandez et al.l 
2010), and it has been known for some time that 
MOND mass-to-light ratios correl ate with color in the 
expected sense in spiral galaxies (ISanders fc Verheiien l 



1998 
2012 ) 



ISanders k McGaughl [2001 : iFamaev fc McGaugl. 



The velocity dispersions given in Tableware genuine a 
priori predictions subject to observational confirmation 
or refutation. Those listed in Table Q] are also genuine 
predictions in the weaker sense often used in the scien- 
tific literature. Though observed velocity dispersions are 
known for the seventeen dwarfs in Table [TJ these do not 
inform the computation of the MOND velocity disper- 
sions in any way Everything follows from the application 
of the equations given in $2] to the photometric data: the 
predicted velocity dispersion follows irrespective of the 
observed one. 



We thank the referee for constructive comments and 
for suggesting the consideration of a narrower range of 
mass-to-light ratios. 
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